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CHAPTER I 


INTRODUCTION 


Business men the world over, have long recognized 
that records and reports of the physical quantities and 
costs of inventories, including their balances, changes and 
relationships, are essential to managerial appraisal of past 
business performance and to planning and control of future 
operations.* The objective of inventory management is to 
maintain the investment in inventories at the lowest amount 
which is sufficient to meet production, sales and financial 
requirements of the enterprise. The inventories must be 
adequate to maintain an efficient level of operations and to 
meet, within reason, the needs of customers. However, it 
must not be greater than is necessary to meet these require- 
ments because of interest costs, the cost of handling and 
storing excessive quantities of inventory, the dangers of 
adverse price changes and obsolescense, and the increased 
exposure to physical detennenationte Thus the requirement 
for an accurate inventory determination is a prerequisite to 
effective inventory management. 

For most business enterprises, the determination of 
inventory assets presents no unusual problems since most 
inventory assets can be physically ascertained. However, 
in certain businesses such as the manufacture of ehemucadls: 


major portions of their inventory assets are in bulk raw and 





finished materials. The problem of accurately determining 
inventory balances of bulk materials has continually plagued 
business managers because such materials are not susceptable 
to recognized mensuration techniques. The quantifying of 
bulk inventories is especially difficult when the economics 
of handling and storing the bulk materials dictate that they 
be stored in stockpiles which have a great degree of surface 


: sponte : 3 
non-conformity and no constraining boundaries. 


Purpose _ and Scope 


The purpose of this thesis is to determine, by experi- 
ment, the minimum measurements necessary to estimate the con- 
tent of irregular solids within a prescribed error uSing nu- 
merical analysis methods. The work encompasses a laboratory 
study utilizing granular ammonium nitrate for construction of 
stockpiles having varying configurations. The study and the 
results of the laboratory tests are to be applied to determine 
the reserves of ammonium nitrate stored in the bulk storage 
warehouses of the Cooperative Farm Chemical Association at 
Lawrence, Kansas. 

In the laboratory experimental project, a digital 
computer program and a manual computation method are de- 
veloped to estimate the content of stockpiles of granular 
Maberzals. Additionally, ae angle of repose of coated 
ammonium nitrate granules is ascertained and compared with 


the accepted value for this characteristic. 


Historical Perspective of Volume Calculations 


For centuries, mathematicians have known the formulas 
for calculating the volume of linear and non-linear solids. 
These formulas require explicit measurements of the length 
of a line or the size of an angle. Examples of such volumes 
are spheres, cones, pyramids and solids of revolution such 
as ellipsoids and cylinders. Volume determinations may also 
be based on implicit measurements. This field encompasses 
the science of integral calculus wherein volumes are computed 
by integration of either a known analytic function or a fin- 
ite numerical difference expression of the area or one linear 
dimension as a function of the other two dimensions. However, 
very little application has been made using the numerical 
difference approach since the amount of. computational work is 
large by this method. Thus, until the advent of the digital 
computer, the best course of action was to estimate the volume 
of irregular solids using the classical explicit volume form- 
ulas and to accept the error inherent with this method. The 
present study is an investigation of the benefits to be 
gained through the use of numerical integration of various 
numerical difference expressions of the area or height of a 


stockpile. 


Problems ,Associated with Stockpile Volume Estimation 


Undoubtably the biggest problem in measuring stock- 
piles is their size. This characteristic coupled with the 


irregular shape of the stockpile surface precludes normal 
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measurement methods based on geometrical formulas. Addition- 
ally, if the stockpile is composed of granular material it as 
Mmipractical, if not ampossibie. vo takercepmen Measurementc 
over the surface because of the fluidity of the granules. 
Consequently, the practice of surveying the surface with a 
imeansit came into vogue in anmeattempt to eStimatomunowaeiccus 
@emthe stockpiles Thissprocedure is quite laborious and tie 
Meetnae, Of the volume calculations leaves something ste) we 
Cieeared . 

Aerial photocontouring is feasible for outdoor stock- 
piles. Results by this procedure are fairly accurate but 
favs method 1s expensive and takes considerable time due to 
the limited number of photocontouring firms. 

Photocontouring on the other hand is impractical for 
femy small outdoor industrial stockpiles. Needless to Say af 
Beenearly impossible to photocontour indoor stockpiles due fo 
factors such as size of the warehouse and environmental con- 
ditions such as dust. 

Another problem associated with determining the con- 
tent of stockpiles is the amount of time and cost involved 
in taking the measurements and performing the necessarily cal- 
culations to fix the inventory quantity. In exercising sound 
inventory management control, engineers and businessmen 
Seould be xnowledgeable of the alternatives available in terms 
of inventory accuracy versus the time and costs required to 


merain that accuracy. 








CHAPTER aan 


NUMERICAL ANALYSIS METHODS 


Numerical Integration 

Integration may be defined as a process of summation. 
Numerical integration is defined by J. B. Scarborough "as 
the process of computing the value of a definite integral 


5 This 


from a set of numerical values of the integrand". 
process iS sometimes called mechanical quadrature when 
applied to the integration of a function of a single vari- 
able; when applied to the calculation of a double integral 
of a function of two independent variables it is called 
mechanical cubature. 


For numerical integration, the problem is to estimate 


the numerical value of the integral: 


b 
I(x) = | f(x) dx (2.1) 

a 
when f(x) either is too complicated an analytical function 
to permit the integration in an analytic manner or when the 
value of f(x) is given only in tabulated form. One method 
of solving the problem is to expand f(x) in an infinite 
series, the individual terms of which can be integrated. 


Thus if 


f(x) =) A, ff, (x) ney 
n=0 





= 
a 
a 





and one can determine the constants 
b 
c,=} B, (x) dx 
a 
then 


b caD 
I(x) = | £(x) dx = yi ca (2.3) 
2 n=O 


Such a method is practical, provided this series converges 
mrerdive 1-.6€., Only a few terms need to be calculated toe 
determine I(x) to any desired accuracy. 

The problem of numerical integration is solved for 
the general case by representing the integrand by a suitable 
polynomial, for a given interval and subsequently integrat- 
ing the polynomial between the desired limits. The poly- 
nomials used for this purpose are called interpolating poly- 
nomials. They are derived using the methods of finite dif- 
ferences. Finite difference methods are covered very exten- 
Sively in the literature of numerical analysis. An excellent 
presentation is that by Scarborough. ° An important property 
of the usual power series interpolating polynomial is that 


th order are con- 


the n'” gifferences of a polynomial of the n 
stant when the values of the independent variables are taken 
waearithmetic progression, that is, at equal intervals. This 
property makes it possible to represent any function, includ- 


ing one in tabular form, by a polynomial if its differences 


euesome degree become constant or nearly constant. 








When the tabulation is made at equal intervals, the 
accuracy of the quadrature formula depends upon the size of 
the interval, the limits of integration, and the degree of 


the interpolation formula. 


General Quadrature Formula 

Scarborough and others present a derivation of Newton's 
formula for forward interpolation which is the basis for de- 
veloping a general quadrature formula for equidistant ordi- 
fetes.( °°?” 


This general quadrature formula is: 








(2.4) 
Xotnh “2 3 “2 AY 
} ae ay. ee en ee eo ee 
xX 
0 
na ae Ay 6 7 en 2. 4 Yo 
org Oe ear se gee 3 eee 
AC eo een 2 Ay 0 
7 ee ee 
Usa So Tan 2 d°y 9 
i reer en ee ee = ON) a | 


where n is the number of equidistant intervals of width h, 
and A is the difference operator on the value of the vari- 
apie y. 

The notation used for the difference operator in the 


general quadrature formula can be explained as follows. If 








Yor Yy Yo °° Yn denote a set of values of any function 

eso CHEN Y, = Yom Yo — ¥qr Vs = )¥eemees Yn nt oe 

soltedethe first differences of the variable y. Denoting 

these differences by Ay), Ay,, Ayo, etc., we have i 

4 a PO? Ay, = Yo oa Yq? AY,_4 = vn ~ J eee ye ss Yn = yo 
The differences of the first differences are called the 


second differences and are denoted by ee Lave ebes 


Bg = 4894-77 OY G0 = Ys ear = a 


Z 
Mya >  8Y¥5 =F AY) = ae ey 


In the like manner, the third differences are: 


3 2 ep 
Sig = Yt ~ VE 6 7 Ye > 2)o - =a ae 


5 = 2 2 
AX¥, = AY - OY] 


Va @¥ae= 82 = ae 


Mimerrecaleanaiysis Iirerature 1. freplere With G@taaedce 
and horizontal difference tables. The use of these tables 
depends on which method of interpolation is to be used. [In 
this study, forward interpolation will be used since the 
fmemlated Gata values begin and end with a zero value over 
the range of integration. 

From the general quadrature formula, specific quad- 
rature formulas may be obtained by varying the value of n. 

Three quadrature formulas have found widespread use 


due to their simplicity and accuracy. These formulas are; 








Trapezoidal Rule, Simpson's Rule and Weddile’'s Rule. They 
are obtained from equation (2.4) by letting n = 1, 2, and 6 


respectively. 


Trapezoidal Rule 


Substituting n = 1 in equation (2.4) and neglecting 


all differences above the first, one obtains 


h 
Ya BY y¢ ee Cya oe may nel ee 
: (25) 
The geometric significance of this formula is that a given’ 
function y(x) is replaced by n/1 arcs (straight-line segments), 
i.e., first degree polynomials. It is exact if y = f(x) is a 


polynomial of degree 1. 


Simpson's Rule 
Substituting n = 2 in equation (2.4) and neglecting 


all differences above the second, one obtains 


Xp*+nh rn 
i ydx = 3 | yo i HH, Ynen*2024%4.. naa 
0 


(2.6) 
Simpson’s Rule corresponds to the representation of the given 
function y(x) by n/2 arcs of a second degree polynomial. It 
is exact if y = f(x) is a polynomial of degree 2 or less. 
While it is not readily apparent, Simpson's Rule is exact if 
y = f(x) is a polynomial of degree 3 because the third dif- 


ference ( a*y0) is zero and the first term dropped involves 


4 
F Yo" 
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Weddic's Rule 
Supstitucine mm — 6 in Cqudtions (2. = vaneance ccETl. 
all differences above the sixth, one obtains 


Xotnh 3 

J ydx = TOHhlygt5 (yy Yo Younes | ee see 

Xo ore 

eC eens) YS 52 eo)” Vomen oer 
C227) 

The geometric interpretation of this rule is that a given 
function y(x) is replaced by n/6 arcs of a fifth degree poly- 
nomial. This rule is exact for fifth-degree polynomial or 
lower. This formula requires at least seven consecutive 
mabues of the function, or, stated alternatively, fhe inte— 
gration interval must consist of six subintervals of equal 


miath h.« 


Volume Computations 


To calculate the volume of any solid using quadrature 
formulas, it is only necessary to establish an x-y grid sys- 
tem and to determine the values of the Z coordinate Conte. 
sponding to each point of that system. Figure (2.1) shows 
the graphical representation for this method of mechanical 
cubature. Thus, if a solid is divided into cross sections, 
each ordinate of that cross section may be measured and tabu- 
lated. These values may then be used in any acceptable quad- 
rature formula such as Simpson's or Weddle’s formulas to 


determine the area, A of each cross section. 


cs’ 
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Figure 2.| 


Graphical Representation of 
Mechanical Cubature ona Irregular Solid 





x 
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b * 
Ls J Z dy (255) 


Subsequently, the cross sectional areas, A may be used as 


cs’ 
the ordinates of the length function and themselves be inte- 


grated using any acceptable quadrature formula. 
d 
Volume = J Aw, 4x (259) 


In this manner it is possible to compute the volume of a solid 
nine ways if one is limited to the three most frequently used 
quadrature formulas based on equal spacing of the independent 
variable. Generally, the accuracy of the double quadrature 
method depends on the number of data points and on the shape 
of the surface of the solid. An important characteristic of 
all the quadrature formulas based upon an interpolating poly- 
nomial is that the formulas are numerically exact for those 
instances in which an né? order polynomial y(x) is used and 
for which instance the actual function being approximated is 


of order n or less. 





CHAPIER EE 


EXPERIMENTAL PROCEDURE 


Apparatus 


The experimental apparatus shown in Figure (3.1) was 
utilized to construct and measure the different configura- 
tions of stockpiles. 

The test site consisted of a 4" x 8' sheet of 3/4" 
plywood mounted on a large stable table. One inch ruled co- 
ordinate paper was then attached securely to the surface of 
the test bed. A rectangular area 24" x 60" was ruled on the 
coordinated paper to establish the boundary of the experi- 
mental stockpiles. These boundaries were used as the ref- 
erence lines in measuring the size of the stockpile at its 


\ 


base. 


Conveyor System 


An overhead track system was devised to simulate an 
overhead conveyor system. A small trolley, carrying a steel 
funnel, was allowed controlled movement along the overhead 
track. The distance from the funnel spout to the test bed 
was designed to be eighteen (18) inches. The model conveyor 
duplicated the formation of stockpiles as they appear in the 


plant warehouse of the Cooperative Farm Chemical Association. 


Measuring Device -- Depth Gauge 
Two laboratory point-depth gauges, Lory Type - A, manu- 


factured by Leupolo and Stevens Instruments, Portland, Oregon 


ae 
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Experimental Apparatus for 


Ceomstructing Laboratory Stale Stockpiles: 
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and capable of measuring distances to an accuracy of .001 
foot, were mounted in two "T" frames which were constructed 
from two inch angle iron. This mechanism permitted the 
measurement of the stockpile surface at one inch intervals 
both across and along the length of the stockpile. The 
gauges were leveled and calibrated for each series of cross 


sectional measurements. 


Materials 

The ammonium nitrate used was taken from the standard 
80 lb. bags supplied by the Cooperative Farm Chemical Associ-~ 
ation. A standard 80 1b. bag of ammonium nitrate was de- 
termined by repeated measurements to contain an average of 
2929.0 cubic inches. The initial density was determined to 
be 47.19 lbs./fcu.ft. Frequent checks were made to determine 
if the stockpile compacted during the four weeks duration of 
the experiment. These measurements indicated that possibly 
avery slight compaction occurred mostly near the center re- 
gions of the stockpile. However the change in height of the 
pile at all positions was within the experimental error in- 
herent in the height measurements. Therefore, it was con- 
cluded that the density change was not sufficient to warrant 
its inclusion in the volume calculations. The original volume 
was therefore used as the basis for all calculations during 


the experiment. 


Miscellaneous Equipment 


A stone jar was used to measure the volume of the 


ammonium nitrate used in the investigation. The volume of 
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the jar was determined by measuring the volume of the water 
necessary to fill it. A calibrated 1000 cc cylinder was used 
to fill the jar. The capacity of the stone jar was deter- 


mined to be 842.0 cu.in. 


Measurement Procedure 

A preliminary series of measurements was conducted 
using a single small conical pile of ammonium nitrate granules, 
to test the proposed measuring procedure. Initially, it was 
proposed that the tip of the depth gauge would contact a 
granule on the pile surface and thereby determine the height 
of the pile at that point. Measurements on the small pile 
using this procedure proved difficult and gave inaccurate re- 
sults. This inaccuracy was caused by the irregular arrange- 
ment of the ammonium nitrate pellets or granules. The tip 
would displace a granule and leave the tip in the air above 
the pile. After repeated measuring attempts, a technique was 
developed whereby the sharp needle-like point of the measure- 
ment probe penetrated the stockpile surface 1/16 of an inch 
(approximately one pellet diameter). This method resulted 
in consistent height measurements. 

The experiment began with the filling of the funnel 
with the ammonium nitrate granules. The nitrate was then 
allowed to drop continuously from the funnel conveyor in the 
overhead track system, falling the eighteen inches to the 
test bed. When sufficient granules had been deposited, such 


that the base of each sub-stockpile reached the ruled 





Lay, 


boundaries, the funnel conveyor was then moved twelve inches 
to a new position. This procedure was repeated until the 

80 lbs. (2929 cubic inches) of ammonium nitrate was deposited 
on the test site. 

Pigure (3.2) shows the configuration of the ‘undis- 
turbed"' stockpile. The four concial shaped sub-piles were 
formed to approximate the typical storage situation in the 
warehouse. The conveyor system was then removed from the 
test site and the height measuring mechanism installed at one 
end. After calibrating the height gauge, measurements were 
taken over the surface of the stockpile at one inch intervals 
in compliance with the requirements of the quadrature formulas. 
In addition to the 1525 height measurements, 61 width mecure= | 
ments were taken of the outline which the base of the stock- 
pile made on the ruled paper. The angle of repose was also 
measured at various cross sections around the stockpile us- 
ing a jointed protractor. Table (3.1) tabulates these measure- 
ments. 

Upon completion of the measurements on the "undisturbed" 
stockpile, a second stockpile was created by selectively re- 
moving some material from the above mentioned "undisturbed" 
pile. It was felt that this would represent the configura- 
tion of a "working" or "disturbed" stockpile. A child's toy 
shovel was used to remove small amounts of nitrate from several 
positions at the base of the undisturbed pile. This process 
Simulated quite well the random removal of the ammonium ni- 


trate by a tractor loader. The amount of 842 cubic inches of 





NS 





Beeure 3.2 


Undisturbed Laboratory Scale Stockpile. 
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TABICE eo - 
ANGLE OF REPOSE AT PEAK 


CROSS SECTIONS AS MEASURED BY PROTRACTOR 
Undisturbed Disturbed 
Stockpile Stockpile 
Cross Section |Tert Right L Right 
Side Side i Side 








mverage Angle 


of Repose 
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ammonium nitrate was removed from the base around the entire 
periphery of the original stockpile. Pigure (3.3) shows the 
"disturbed'' stockpile. After the removal of the 842 cubic 
inches of material no attempt was made to alter the shape of 
the pile. As can be seen from Figure (3.4), the surface 
tends to be linear and the peaks have regressed into ridges. 
The "disturbed" stockpile was then measured in the same man- 


ner as was the "undisturbed" stockpile. 


Computational Models 


An explanation of the various computational models 
employed to ascertain the volume of the "undisturbed" and of 
the "disturbed" stockpiles will now be given. These models 
range in complexity from the simplist in which the stockpile 
cross sections are assumed to be isosceles triangles to the 
most complex model wherein the depth measurements on a one 
inch square grid are numerically integrated using Weddle'’s 
Rule. The computer programs for each model were devised to 
allow parameters such as angle of repose, and the number of 
measurements in the x and y directions of the x-y base grid 
to be considered as variables, to be read in as data. This 
flexibility permitted a study of the influence of these 
parameters on the volume calculation. 

The computer programs for the computational models 


used in this study may be found in Appendix A. 


Model # 1 
A computer program written in Fortran IV language 


was developed to compute the volume of stockpiles assuming 
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Disturbed Laboratory Scale Stockpile. 
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Figure 3.4 


Disturbed Laboratory Scale Stockpile 


Depicting the Change in the Angle of Repose. 
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isosceles triangular cross sections. The program calculates 
the cross sectional areas utilizing the angle of repose, a, 
as a parameter. Additionally, in computing the volume by 
numerically integrating the cross sectional areas over the 
length of the stockpile, three different quadrature formulas 
were used. The number of cross sections was also varied to 
observe its effect on the volume calculations. 

The equation for the area of the isosceles triangular 
cross section is shown below for the circumstance in which 


the angle of repose, a, and the base b are known 





tan a = h/s 
Cross Section Area = 1/2 bh 
= (1/2)(b)(B tan a) 


2 
b : 
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Model # 2 | 

A computer program was developed to estimate the volume 
of a stockpile by numerical integration in two aiiagetionen 
The program calculates three different areas at each cross 


section utilizing the three quadrature rules: 


1) Trapezoidal Rule, 
2) Simpson's Rule, 


3) Weddle'’s Rule. 


Also the number of depth measurements used in each cross sec- 
tional area can be designated externally and used as a model 
parameter. The program calculates the volume by numerical 
integration of the cross sectional areas previously obtained. 
For this phase of computation the number of cross sections 
and particular quadrature rule to be used can be selected 


externally and used as model parameters. 


Model # 3 

This computer program is a modification to Model #1 
and was written to specifically compute the volume of a stock- 
pile which is banked against a retaining wall. 

The cross sections are assumed to be isosceles triangles 
but with a portion of the triangle removed. The net area of 
each cross section is determined by computing the area of the 
isosceles triangle and then subtracting that part of the isos- 
celes triangle which projects beyond the retaining wall. 

Using the angle of repose as a parameter, the program 


proceeds to compute the volume of the stockpile by numerical 





integration of the previously obtained cross sectional areas 


Pome gene three quadrature rules, 


The equation for the area of a typical cross section 


is derived below for a banked stockpile. 


h 
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CHAPTERS IY 
ANALYSIS OF EXPERIMENTAL RESULTS 


General 

The results of the experimental study outlined in 
Chapter III will be presented in three parts. Part Ais 
concerned with the computation of the volume and the deter- 
mination of ane apparent angle of repose for the "undisturbed" 
stockpile. Part B describes the same results but for the 
"disturbed" stockpile. Part C describes the results from the 
determination of the true angle of repose for coated ammonium 


Haetrate granules. 


A. Undisturbed Stockpile Results 


Computational Models #1 and #2 were used in obtaining 
the volume and apparent angle of repose approximations. 
Model #1 

Table (4.1) tabulates the per cent error in the cal- 
culated volume obtained from computational Model #1 using 
aS parameters; 

1. number of isosceles triangle cross sections 

2. quadrature rule 

3. apparent angle of repose. 

Because of the mathematical form of Model #1, the apparent 
angle of repose, a, is an independent parameter. The volume 
of the stockpile computed by the model is dependent upon the 


value selected for a. If the volume of the ammonium nitrate 
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TABLE. 421 


PER CENT ERROR SNe CALC UAT ED? VOLUME “POR 
UNDISTURBED STOCKPILE USING COMPUTATION MODEL #1 







Per Cent Error in Calculated Volume 
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as measured in the volumetric jar is taken as correct, then 
the per cent error in volume estimated by Model #1 can be in- 
dicated as a function of a. The results are depicted in 
Figure (4.1). 

Before analyzing these results it would seem appro- 
priate again to describe the manner in which the measure- 
ments were taken. A rectangular grid boundary, 24 inches by 
60 inches was drawn around the base of the stockpile. Thus, 
Sixty-one cross sections could be selected, at a one inch 
spacing, the first and last cross sections having zero cross 
sectional area. By measuring the distance between the base 
of the stockpile and the side boundary for each cross section, 
the base length of each individual cross section was deter- 
mined. The accuracy of these measurements was approximately 
1/10 of an inch because the scale was ruled in tenths and 
because individual particles were not much smaller than 1/10 
inch. 

From Figure (4.1a), it may be observed that the three 
quadrature rules give essentially the same error curve for 
volume determinations when the maximum number (61 ) of sec= 
tional areas is used in the computation. From the graph, 
the average apparent angle of repose for coated ammonium ni- 
trate granules was determined to be approximately 29.7 de- 
grees. The angle was taken to be the point at which the per 
cent error in the volume determination is zero. This appa- 
rent angle of repose would be used in Model #1 to estimate 


the volume of an actual undisturbed stockpile having 
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Error in Volume Estimation by Model Number | 
for an Undisturbed Stockpile 
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essentially the shape of the one used in this particular 
test. The stockpile could consist of one or more peaks all 
essentially in line. 

Further study of Pigure (4.1) reveals that the error 
curves for the three integration rules used in Model #1 di- 
verge as the number of cross sections is decreased. [In order 
to study the effect of varying the number of cross section 
measurements on the volume calculation, additional volume cal- 
culations were made for 21, 16, 11, 6, 5, 4, and 3 cross sec- 
tions. These results together with the results for 61, 3l, 
13, and 7 cross sections from Table (4.1) are listed in 
Table (4.2). 

Table (4.2) presents the per cent error in the calcu- 
lated volume utilizing the apparent angle of repose of 29.7 
degrees. The results are presented graphically in Figure (4.2). 
The per cent error in the volume behaves as a variable under-= 
going damped oscillation. The error is large in magnitude 
for a small number of cross sections becoming damped to a 
small constant error as the number of cross sections ap- 
proaches one-half of the total available. Specifically, the 
figure indicates that the error curves converge at approxi- 
mately 31 cross sections and are quite close when 15 to 31 
cross sections are used. Thus, to compute the volume of a 
stockpile having essentially the same shape as the stockpile 
used in the present laboratory study one should use at least 


20 cross sections and possibly as many as 31 if the volume is 
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TABLE 4.2 
PER CENT ERROR IN “CALCULATED VOLUME FOR AN 


NEPARENT ANGLE OP REPOseeeGe7o.7 DEGREES POR 
MODEL #1 


Per Cent Error in Volume 
Cross Section/ Trapezoidal Simpson's Weddle's 
Interval Size* Rule Rule Rule 


0.00 
-0.52 


o/ 1 
S/2 
21/3 
16/4 
fs: / 5 
11/6 
@/ 10 
oy 12 
By/ 1.5 
4/20 
3730 





* The notation 61/1, 31/2, etc., indicates 1 
interval containing 61 points, 2 intervals 
containing 31 points each, etc. Note that 
the terminal point on one interval is the 
Thitial point on the nexterarceryal, thus 
(points -1) x no. of iantervamee— 960 for all 
Se ts.. 
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Figure 4.2 


Error in Volume Estimation by Model Number | 
| for an Undisturbed Stockpile 
(Angle of Repose = 29.7°) 
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to be estimated by Model #1. The appropriate apparent angle 
of repose would, of course, be 29.70 degrees. 

It is believed that the results of the above labora- 
tory model can be scaled to the dimensions of commercial 
stockpiles by considering the number of peaks in the commer- 
cial undisturbed stockpile. Thus, assuming the minimum num- 
ber of cross sections is 20 for the 4 peaked pile studied 
above, (see Pigure 3.2) one might generalize the results to 
a rule that a minimum number of 5 cross sections per peak, 
i.e., 20 sections/4 peaks, should be used to estimate the 
volumes of undisturbed multi-peaked stockpiles of the type 
studied herein. 

It is to be noted that the Trapezoidal Rule gives a 
relatively low error level over the greatest range of cross 
sections. The Trapezoid Rule calculates the volume with an 
accuracy of = 1.0% for 13 or more cross sections while 
Simpson’s Rule calculates the volume to this Same accuracy 
for 21 or more cross sections. Weddle’s Rule converges 
Slower and therefore requires a greater number of cross sec- 
tions to compute the volume within u 1.0%. 

Model #2 

Table (4.3) tabulates the per cent error in the cal- 
culated volume obtained from computational Model #2 which 
involves two-way numerical integration. The following para- 
meters were varied in this study: 


1. number of cross sections 
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PERCENT ERROR IN CALCULATED VOLUME POR UNDISTURSED 
STOCKPILE USING COMPUTATION MODEL #2 


ee Number of Depth Measurements 
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2. number of depth measurements per cross section 

3. quadrature rule for cross sections 

4. quadrature rule for length 

Because of the difficulty of conveniently represent- 
ing graphically all of the data containedin Table (4.3) only 
certain extreme examples are plotted herein. These corre- 
spond to the three curves which give consistently the small- 
est error in the volume calculation when uSing both the maxi- 
mum number of points (25 points) in the cross sectional area 
determination and the minimum number of cross section points 
(3 points). These results are plotted as Figure (4.3). 

It may be seen that the error curves have the same 
damped oscillation feature observed in Model #1 with the 
error becoming relatively constant when approximately one 
half of the total possible cross sections on a one inch 
Spacing are used, regardless of the number of depth measure- 
ments per cross section. AS might be expected, the error is 
greater for the smaller number of depth measurements used per 
cross section. 

The results of Figure (4.3) show that a remarkably 
low error is incurred by using as few as three points per 
cross section when a large number of lengthwise cross sec- 
tions are taken. This error, less than 3%, is not appreci- 
ably greater than that obtained with a large number of cross 
sections and the assumption of an isosceles triangular cross 


section using Model #1. 
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Figure 4.3 


Error in Volume Estimation by Model Number 2 
for an Undisturbed Stockpile 
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This observation seems to indicate that, for Model #2, 
one need not use a large number of data points in both the 
width and length dimensions on an undisturbed stockpile to 
obtain a satisfactory estimate of stockpile volume. Instead, 
one is probably justified in using a rather elementary ap- 
proach in obtaining cross sectional areas provided a sufficient 
number of these cross sections is used for the lengthwise in- 
tegration. 

An analysis of the calculated areas for two typical sec- 
tions was made to determine the best method of obtaining 
the area of any cross section. Table (4.4) tabulates the 
results of the area determinations for two typical cross sec=- 
tions for the three quadrature rules. Using the maximum 
number of depth measurements, all quadrature rules computed 
the areas for the typical cross sections within .50%. 

Using the results obtained by applying Weddle's Rule 
to the maximum number of cross sectional depths measurements 
as the "correct"! answer and then comparing the results ob- 
tained from both Simpson's and the Trapezoidal quadrature 
Rules based upon only three depth measurmments, it may be 
observed that the Trapezoid Rule gives the lowest error. 

In a study similar to that shown in Table (4.4), but 
involving more cross sections, it was revealed that the 
errors in area estimation by the Trapezoidal Rule were oppo- 
Site in sign for cross sections of high and of low area. 

These errors tended to cancel when used in the lengthwise 


integration formulas, thus further substantiating the use 
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TABLE 4.4 


COMPARISON OF COMPUTED CROSS SECTIONAL AREAS 
FOR UNDISTURBED STOCKPILE 


Cross Sectional Area Coe 


Number of Depth Quadrature Number 13 Number 7 
Measurements per 
Cross Section Rule (large) (small) 


25 Weddle's 74.66 D7 

25 Simpson's 74.62 ey . 98 

25 Trapezoidal 74.62 27 .86 
ren 

3 

3 








a aE | a 
Trapezoidal err OO Sloe 
Simpson's 96.00 41.86 


Considering the 25=-point Weddle's Rule 
values to be the "correct" areas for 
each cross section the error at selected 
cross sections are as follows: 





Cross Section #13 


: a (72.00 — 74.66) 
Per Cent Error in Area 74.66 LOOR® 


-3.56% 


. or Error in square inches = -2.66. 


Cross Section #7 


(ieee 26 97) 


Per Cent Error in Area = 4707 


L0O0.0 = +1232 


or Error in square inches = +3.42. 
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of the simple Trapezoidal Rule for obtaining cross sectional 
areas. 

Table (4.5) summarizes for Model #2 the best alton 
tives in keeping with the basic purpose of this study - namely 
to determine the minimum number of measurement points and 
their location in order to best approximate the volume of a 
stockpile within the prescribed error range. The table pre- 
sents the alternatives for which the volume determination on 
the undisturbed stockpile, computed by Model #2, is less than 
= 3%. 

From the observations of the volume calculation methods 
for the undisturbed stockpile it may be stated that both 
Models #1 and #2 give excellent results when utilizing a 
maximum number of data measurements. Also, under appropriate 
circumstances, a much smaller number of measurements can be 
used in the models and still maintain the error levels within 
a tolerable range. Specifically, one is justified in using a 
rather simple geometrical representation in the cross-wise or 
cross sectional area determinations for elongated undisturbed 
stockpiles. It is possible to determine the cross sectional 
areas with a high degree of accuracy from as few as 3 to 5 
measurements, with approximately 5 cross sections being taken 
per peak over the length of the stockpile. 

If this procedure is followed using the Trapezoid or 
Simpson's quadrature Rule, then the stockpile volume should 
not be in error by more than = 3%. Further refinements ager 


of course, possible using additional measurements but 
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TABLE 4s 


ALTERNATIVE MEASUREMENTS ON AN UNDISTURBED STOCKPILE 
FOR VOLUME ERROR RANGE OF - 3.0% 


Number of Depth Number Total Number Double 
Measurements per of ! of Quadrature 
Cross Section Cross Sections Measurements Rule* 


WDSN 
WDTR 
WDWD 





Note: 


Pigures in pasememeetsmane actual 
number of measurements required 
Since the end ©f ourer boundary 
measurements are always zero. 


*The Code for the Double Quadrature Rule is as follows: 


WD Weddle's Rule 
SN = Simpson's Rule 
TR = Trapezoidal Rule 
thus WDSN means Weddile's Rule used in integrating 


the cross section and Simpson's Rule used to inte- 
Brate sever“Ehe Length. 
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accuracies greater than 1% are not to be expected. On the 
other hand, a very rapid approximation method using only 3 
points (in reality only 1 point) per cross section and 3 


cross sections per peak will give an error of less than 10%. 


Bea Disturbed Stockpile Results 


Computational Models #1 and #2 were again used in 
obtaining the volume and apparent angle of repose of the 
disturbed stockpile. The method of analysis was essentially 
the same as that for Part A of this chapter. 

Model #1 

Table (4.6) presents the per cent error in the cal- 
culated volume obtained from computational Model #1 using 
the following as parameters; 

1. number of isosceles triangle cross sections 

2. Quadrature rule 

3. apparent angle of repose. 

These results are shown graphically by Figure (4.4). It can 
again be observed that the error in the volume estimation by 
the three quadrature rules is essentially the same when the 
maximum number of cross sections is used in the computation. 
The graphs show that the apparent angle of repose for a dis- 
turbed stockpile is approximately 25.9 degrees. This value 
closely corresponds to the value of 26.5 degrees used by ene 
Cooperative Farm Chemical Association in their calculations 

on disturbed stockpiles. Further observations of Figure (4.4) 


shows that the error curves are less widely separated for the 





42 


TABLE 436 


PER CENT ERROR IN CALCULATED VOLUME FOR 
DISTURBED STOCKPILE USING COMPUTATION 
MODEL #1 
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range of 7 to 61 cross sections than the corresponding curves 
applicable to the undisturbed stockpile. This implies a 

more nearly constant apparent angle of repose as the number 
of cross sections is decreased from the maximum of 61. This 
may be attributed to the increased homogeneity of the cross 
sections caused by the removal of material from the base. 

By reference to the disturbed stockpile, Figure (3.3) one 

can see that the peaks of the undisturbed stockpile have 
formed into ridges as a consequence of the removal of material 
from around the bottom of the stockpile. This action created 
a relatively planar stockpile surface in the region of the 
ridges. The forming of the ridges from each peak is there- 
fore, visual evidence of the homogeneous structure of dis- 
turbed stockpiles. 

Further volume calculations were made for the same 
additional cross sections as in Part A and the results are 
listed in Table (4.7) together with the applicable results 
of Table (4.6). Table (4.7) presents the per cent error in 
the calculated volume utilizing the apparent angle of repose 
of 25.9 degrees. 

Figure (4.5) depicts the error in the calculated 
volume as a function of the number of cross sections. For 
both the disturbed and the undisturbed stockpiles, approxi- 
mately 20 to 30 cross sections seem sufficient for an accu- 
rate volume determination with possibly slightly more cross 
sections being required for disturbed stockpiles since the 


curve for the disturbed stockpile does not seem to damp quite 





TABLE. 


PER CENT ERROR IN CALCULATED VOLUBS FOR 
APPARENT ANGLE OF REPOSE OF 25.9 DEGREES 
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Figure 4.5 


Error in Volume Estimation by Model Number |! 
for ao Disturbed Stockpile 
(Angle of Repose = 259°) 
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as rapidly as for the undisturbed stockpile. It is not known 
however, how general is this difference in damping between 
the two types of stockpiles. Por Model #1, an accuracy of 

= 3% is attainable with as few as 20 cross sectional measure- 
ments. 

In general, it may be concluded that the disturbed 
stockpiles require a slightly larger number of cross section 
intervals than do undisturbed stockpiles to compute the 
volume with the same accuracy. 

Model #2 

Table (4.8) presents the per cent error in the cal- 
culated volume obtained from computation Model #2 using the 
same parameters for Model #2 as used in Part A. A compari- 
son of these results with those listed in Table (4.3) for 
the undisturbed stockpile shows a general overall increase 
in the error of the computed volume for the disturbed stock- 
pile. This is especially noticeable when the number of 
depth measurements utilized in the computation, is decreased 
to a minimum. The error involving Weddle's Rule using three 
depth measurements iS an exception; the error value is mean- 
ingless with this small number of points since Weddle'’s Rule 
requires 7 points. 

Pigure (4.6) depicts the per cent error in the cal- 
culated volume obtained from the double quadrature rules 
which gave the minimum error calculation. Two sets of error 
curves are plotted; one utilizing the maximum number of 


depth measurements per cross section (23 points) and the 
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TABLES 


PER CENT ERROR IN CALCULATERSV@LUME FOR DISTURBED 
STOCKPILE USING COMPUTATION MODEL #2 


Depth Measurements 


23 Pacman | 7 PTS 
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Figure 4.6 
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other using the minimum number of depth measurements per 
cross section (3 points). The dampning characteristic for 
both sets of error curves is somewhat diminished compared to 
the corresponding curves (see Figure 4.3) for the undisturbed 
stockpile. The most noticeable trait of this plot is the 
relatively large difference between these two sets of error 
curves. The set of error curves using the maximum number of 
depth measurements per cross section shows a = 2% error for 
the range of 31 to 61 cross sections and an overall error 
level of is 6% for 7 to 61 cross sections. On the other hand, 
the set of error curves utilizing 3 depth measurements give a 
minimum error of approximately 10% over the complete range 

of cross sections. 

These results indicate that two-way numerical inte- 
gration over disturbed stockpile is not very accurate using 
a minimum number of depth measurements per cross section un- 
less one uses the Trapezoid Rule to determine the cross sec= 
tional areas and then integrate the resultant areas by either 
the Trapezoid or Simpson’s Rule to obtain the volume. 

It can therefore be concluded that for the disturbed 
stockpiles, which have relatively straight surfaces, the 
Trapezoidal quadrature Rule should be used for integration 
when using the minimum number of depth measurements per 
cross section. The volume error level which results using 
this method is reasonably stable but is of such a magnitude 
(10%) that another method for calculating the volume is to 


be preferred. 
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Table (4.9) presents for Model #2 for a disturbed 
stockpile, the best alternative measurements required to cal- 
culate the error in the volume to an accuracy of + 3.0%. A 
comparison of this table with the corresponding Table (4.5) 
in Part A shows that a greater number of depth measurements 
would be required on a disturbed stockpile to achieve the 
same 3.% error level. 

This study has shown that numerical integration when 
using a minimum number of cross sectional depth measuréments 
resulted in a error level that was 3 times greater than that 
determined by Model #1 using isosceles triangular cross sec- 
tions. It is therefore concluded that the use of the single 
isosceles triangle cross section is justified in approximat- 
ing the volume of disturbed stockpiles. 

In summary it may be stated that Model #1 seems pre- 


ferable for computing the volume of disturbed stockpiles. 


C. Angle of Repose Results 


The angle of repose of a granular stockpile is defined 
as the angle which a plane, tangent to the surface of the 
stockpile, makes with a horizontal plane. In this experiment, 
the undisturbed stockpiles consisted of sections that were 
conical and for this geometrical shape the angle of repose 
corresponds ee the angle formed by the horizontal plane and 
by a line lying in the conical surface and passing through 
the apex of the cone. 

Since height measurements were taken at one inch in- 


tervals, in general, a cross section would not pass precisely 
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TABEEV 427 


ALTERNATIVE MEASUREMENTS ON A DISTURBED STOCKPILE 
FOR A VOLUME ERROR RANGE OF = 3.0% 





Number of Depth Total Number Double 
Measurements per of Quadrature 
Cross Section Measurements Rule* | 
TRUER 
23 61 1403 SNTR 
iN 
SNSN 
23 31 TiS | TRIR 
TRSN 
| dx SiN 
23 7 161 ONTR 
| SNWD 


TRSN 

7(5) 31 317(155) TRTR 
TRWD 

7(5) 13 91 G65 esl TRTR 
SNTR 

7(5) 7 49(35) =a 


Note: 


FigureS in parenbWesis are actual 
number of measurements required. 
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through the apex of any of the four conically shaped peaks 
of the undisturbed stockpile. Therefore, in order to esti- 
mate the true angle of repose it was necessary to interpo- 
late between the cross sections bracketing the cross section 
which passed through the apex. 

A straight line was passed through the depth measure- 
ments, by the method of least squares, for each side of all 
cross sections. The slope of each line was determined and 
the angle less than 90 degrees corresponding to this slope 
was computed. These angles were then plotted versus the 
cross section position. The true angle of repose was then 
taken to be the interpolated angle corresponding to the 
known peaks cross section position. FPigure (4.7) shows a 
typical plot of angles computed in the manner described above 
and corresponding to cross sections near as undisturbed stock=- 
pile peak. The graph shows that the angle of repose is about 
31.6 degrees at cross section location 13.2. Table (3.1) 
lists the observed angle of repose by protractor measurement 
for cross section number 13 as 31.0 degrees. Figure (4.8) 
shows a plot of the height measurements for cross section 13 
together with the straight line obtained by the method of 
least squares. 

The angle of repose was also determined in the same 
manner, for the disturbed stockpile. Figure (4.9) is the 
plot of the angles of repose for the cross sections forming 
the first ridge section (observed at the left in Figure 3.3) 


Since the apex of the conical shape degenerated to a ridge, 





_°? SE 





54 


Figure 4.7 
Angle of Repose for Undisturbed Stockpile 
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Figure 4.9 


Angle of Repose for Disturbed Stockpile 
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it can be seen that it is now possible to have adjacent angles 
of repose varying from 28.0 to 30.2 degrees for a disturbed 
stockpile, the average being approximately 29.5 degrees. 

One may think of the 31.6 degrees angle of repose for 
an undisturbed stockpile as the initial condition of a stock- 
pile. With the passage of time, factors such as weight of 
the granules, vibration, moisture, and size of granules cause 
the stockpile to compact and thereby reduce the angle of re- 
pose. For this experiment lasting approximately four weeks, 


the angle of repose decreased 2.1 degrees. 





CHAPTER V 
APPLICATION OF RESULTS 


General 

This chapter discusses the principles and methods of 
volume estimation for commercial scale stockpiles based on 
the methods developed in the previous chapters for laboratory 
scale stockpiles. The use of a digital computer program and 
a manual computation method utilizing the stockpile base 
measurements and the apparent angle of repose is explained. 
The computer program, Model #1, given in Appendix A, is 
adopted as the most appropriate mathematical model to be 
used in estimating the volume of stockpiles. This mathemat- 
ical model, which assumes that the typical stockpile cross 
section is an isosceles triangle, is the basis for the de- 
velopment of the manual computation forms shown in Appendix 
B. Appendix B contains one blank manual computation form 


and one completed sample calculation. 


Reserves Estimation Using a Digital Computer 


In using Model #1 computer program to compute the 
volumes of commercial scale stockpiles, it is first necessary 
to establish a boundary system in the stockpile warehouse. 

The boundary system can be either painted on the floor or 
associated with the stanchions of the warehouse. The side 
boundaries should be located as close to the base of the stock- 


pile as is physically possible. Also, the end boundaries 
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should be selected after proper recognition is taken of the 
number of cross sections to be used per peak. The total num- 
ber of cross sections taken should be an odd number. Further, 
although not absolutely necessary, it would be wWisePrGrse acim: 
the interval between cross sections so that the cross sec- 
tions all fall on marked scale divisions of the base grid. 
These criteria allow a maximum flexibility in the use of the 
Trapezoid and Simpson's integration formulas, and maximum 
Simplicity in making the measurements. 

In determining the number of cross section base measure- 
ments to be taken, a decision needs to be made as to whether 
the stockpile is classified as "undisturbed" or "disturbed". 
If the stockpile is undisturbed, the general rule is to take 
at least three base measurements per each peak on the stock- 
pile. The resultant error in the estimated volume will be 
within approximately =< 3%. If more than 3 measurements per 
peak are taken the resultant error in the estimated volume 
can be as low as ~ 1%. 

On the other hand, if the stockpile is disturbed, the 
eemeral rule is to take at least 5 base measurements per each 
peak. The error in the estimated volume will be approximately 
in the range of = 3%. If more than 5 measurements per peak 
are taken the error in the estimated volume will approximate 
= 1%. 

Before a final selection is made as to the number of 
measurements to be taken, consideration should be given to 


the rough estimate of the amount of material in the stockpile. 
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Good inventory management practice dictates that the larger 
the stockpile, the more accurate should be the volume deter- 
mination. -Thus, the maximum number of base measurements 
should be taken consistent with the labor costs necessary to 
take the measurements. 

Once the number of base measurements to be taken has 
been determined, the next step is to begin measuring the 
distance between the base of the stockpile and the side 
boundary line. The measurements should be to an accuracy of 
1/10 feet or the nearest inch. The cross section base measure- 
ments are required to be equally spaced over the length of 
the boundary grid and should of course be taken on both sides 
of the stockpile. The measurements should be recorded on a 
data sheet similar to that shown in Figure 1, Appendix B. 

Upon completion of the base measurements, the volume 
and tonnage estimates for the stockpile can be computed with 
a digital computer. The measurement data should be keypunched 
following the instructions given in Appendix A for the Model 
micomputer program. The computer program also requires the 
length and width of the boundary grid together with the appa- 
rent angle of repose to be read in as data. The apparent 
angle of repose must be determined for different kinds of 
stockpile material in the manner outlined in Chapter III. The 
proper formats for all data cards are prescribed in the com- 
puter program. | 

The computer program computes the estimated volume in 


the units selected for the input data. 
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The weight of material in the stockpile is then cal- 


culated using the average bulk density of the material. 


Reserves Estimation by Manual Calculation 


Appendix B contains both a blank manual computation 
form, Figure 1, and a completed manual computational form, 
Figure 2. This form was designed to use the same measure- 
ments as described for the computer program. The form is 
designed in three basic parts. At the top of the form is a 
block for constants. These are self explanatory and are 
used in the remaining two parts for computing the cross sec- 
tional areas, volume, and weight of stockpile material. The 
center section of the form calculates the areas of the cross 
sections utilizing the base measurement distances in Columns 
2 and 3 and assuming that a typical cross section is an isos- 
celes triangle. Columns 7, 8, and 9 of the form are totaled 
and these results form the basis for the volume and tonnage 
determination of the stockpile. 

In the computation form, provision is made to compute 
the reserves by two different means; one using the Trapezoid 
Rule, the other using Simpson's Rule for the lengthwise in- 
tegration of cross section areas. In general, these will 
give close but not identical answers as in the present case. 
It 1s suggested that either the two answers be averaged, or 
the most conservative estimate be selected and used for in- 


ventory management purposes. 





CHAPrER ey I 
SUMMARY AND RECOMMENDATIONS 


in thas thesis an at tenptieneospcensmade to improve 
the volume estimation of commercial stockpiles through the 
judicious use of mensuration and numerical integration 
methods on laboratory scale stockpiles. 

During this study, the basic principle of mensura-~ 
tion was kept foremost in mind, namely, to determine the 
area of surfaces or the volume of solids from certain 
feeriple data on lines and angles: 

Using numerical analysis methods, it was found to be 
advantageous to unite the two computational methods. The 
experimental results show that the areas of all cross sec- 
tions of a stockpile can be readily and accurately approxi- 
mated by an isosceles triangle. Moreover, the area of an 
isosceles triangle can be determined by two simple measure- 
ments, its base length and angle of inclination. The appli- 
cation of either the Trapezoid or Simpson Rule of numerical 
integration to the approximate cross sectional areas re- 
sults in the most accurate volume estimates. The choice 
of the rule is somewhat arbitrary but depends essentially 
on the degree of regularity of the stock pile surface. It 
is recommended that the Trapezoid Rule be used for disturbed 
surfaces and the Simpson Rule for undisturbed surfaces, 
mee., those formed by pouring material directly above the 


pile. Thus, advantage can be taken of both mensuration and 
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numerical integration concepts to satisfy the requirement for 
a minimum number of measurements to attain a relative high 
Becumacy for thesvolume determinacaon. 

The study has shown that, for a laboratory scale 
stockpile, this method of ea ete the volume is realis- 
tic. It is believed that the results of the present Labora- 
tory measurements can be adapted to the dimensions of com- 
mercial stockpiles by considering the number of cross section 
measurements taken per stockpile peak. For practical use, 5 
cross sections per peak are sufficient although 3 cross sec- 
tions per peak generally give errors of less than 10%. 

Finally, if this technique has awakened engineers and 
businessmen to the prospects of making a better determination 
of materials stored in a stockpile, for inventory management 


purposes, the main purpose of this thesis has been accomplished. 
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APPENDIX 


COMPUTER PROGRAM TO DETERMINE THE APPRC..IMATE VOLUME OF A 
STOCKPILE ASSUMING THE STOCKPILE CROSS SECTIONS ARE ISOSCELES 
BR LANGLES. 


Purpose 


The purpose of this program 2eeto determine the aperess 
mmate volume of a stockpile by numerical integration assuming 
meosceles triangle cross sections. A secondary purpose is To 
determine the angle of repose of the material forming the 


stockpile. 


LANGUAGE 


Fortran IV (CIBM 7040 Computer) 
exmMoeOLIC DICTIONARY 


VARTABLE 7 a Lj DESCRIPTION 


XR A I Distance the base of the 
pile is from right bound- 
Beye te. . 

XL A i Distance the base of the 
pile is from left bound- 
age tane:. 


M S I/O Number of measurement 
Points. 
1B S I/O Distance between the end 


boundary lines. 


W S I/O Distance between the side 
BOumGary lines. 


*S - Single variable; A ~ Array of variables 


**XT = Input; O - Output 





Pole ri 


mip PHAR 


N 


VACT 


PpELTAL 


BASE 


AREA 


SUMA 


PCALTR 


PERRTR 


moaAREA 


EVAREA 


VCALSN 


PERRSN 


I/O 
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Inclination angle of the 
stcckpile -— same as appar- 
ent angle of repose. 


Alpha expressed in radians. 


Number of measurement 
Sue. 10ns. 


Acétual volume of the 
SEOcKDa Lee 


Meneth anterval sages 


Width of stockpile at the 
base. 


brea Of -a-Cross Section. 


A holding variable used 
Peomaeccumiuitate the areas 
Grecross Sections for cal= 
Gelating the volume Dy tie 
meapezeraal Ruler 


The calculated volume of 
Pnevetockorile utilizing ime 
fiuapezoidal Rule. 


Per cent error of the cal- 
euleateda volwmpe using tne 
Peewee 2Oraa 1. Rute. 


Pehorries variable used ee 
accumulate the areas of the 
SPQGenmunbered SEOCKDiLe 

cross sections for calculat- 
ang the stockpile volume us-= 
Piss oimpSen’s, Rule: 


Peemolaane Variable used ie 
accumulate the areas of the 
even numbered cross sections 
for caiculating the stock- 
pire volume using Simpsconme 
Rue. 


Pree Ccaleutlared stock pire. 
volume using Simpson’s Rule. 


Per cent error of the cal- 
culated volume using 
Simpson’ ssRuie - 
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SUMTWA S a A hol@inge variabic useadaro 
accumulate the weighted 
areas of cross sections for 
calculating the stockpmle 
volume using Weddle" ss kate: 


VCALWD S O Tne scatculared St o¢ipaee 


volume uSing Weddle s hate 

PERRWD S O Per cent ¢rror of thes 4at- 
Cutated volume uSing Weda@ies 
mie, 


PROGRAM ROUTINE 

This program utilizes the data points (representing 
the distances to the side boundary lines) for each cross 
section. The area of each cross section is then computed 
uSing the angle of repose as a parameter. These areas are 
then integrated by the Trapezoid, Simpson's and Weddle's 
meadrature rules for a set interval Sizé to obtain the cal- 
culated volume of the stockpile. The per cent error in the 
Salculated volume is then computed. The angie of repose is 


determined when the per cent error iS zero. 








MOM 
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THIS PROGRAM DETERMINES THE APPROXIMATE VOLUME OF A STOCKPIL®E 
ASGUMING THAT THE CROSS SECTIONS S8— 1SOCELES Tir eee. 


BIWeENSION XL(100) ss XkR(1007"s AREATIC?C) 
READ (555) M 

MORMAT ( h2) 

WRITE ( 656) 

mMORMAT (15Xs 20HM = Niger DATA PISse//)} 
WRITE (657) M 

Pome! (19Xs 127/70 

MEAD (5520) (XLII) » Ameri) s91 = le M) 
moRMAT (2F 1Oe2) 

Vetwe (6s 1 ) 

FORMAT (Xs 17APNEUT DA eorINTS 

Ret Lent .6 61.3.) 

FORMAT (15Xs 3HNOes 15X9 2HXLe 15X9 2HXR//) 
WRITE (6212) (¢( Iso XLC1)e XR{(I1) =» 1 = le M) 
FORMAT ( 15Xe [259 1°Xs Fide2s Bis FIGeZz// } 
READOW 5610) Ls Wo AIRE tee ov CI 

FORMAT (2F 10 «26 FGmeee pom cs ob |e 

WRITE (6914) 

FORMAT (15%. PORTNOY BA IAREONSCTANTS//) 
Mikel | Ct Gredee } 

FORMAT(15X%> lH SOX s TriWer Ove SHALPHAs 1O0Xs LHNs 90K eee 
WRITE (6597) LP WRAL PHA sth suey ee] 
MORMET(I3Xs F5e2s CXs Pome se OK oe 2s CXS oe eo F1lQe2//) 


ASSUME CROSS SECTICNS ARE ISOCELES TKIANGLES FOR COMPO | Gea 


REAL Lo»N 
ALPHAR = 3414159 * ALPHA / 180.0 
C = (SIN(ALPHAR}/COS(ALPHAR) )/4¢0 


DELTAL = L/(N-1.)0) 

Wied (601 6) 

FORMAT (15X+ 17HINTERVAL = DELTAL//) 
(fete AOsl9) DELTA 

Boel) (22Xs FouZzyye 


J=zn 

Peewee, = 1s 2 

BASE = W —- XL(1I) -— XR(1) 
AREAL] ) = CC *htSr Gee 


WRITE (69500) 

FORMAT (25Xs- 22HAREA OF S@- Ce PCT IONS//} 
WRITE (63600) (Ie AREA( 1)» I = le M) 
PORPeAr«15xX= l25 Fille) 


STOCKPILE VOLUME COMPUTED yi THE TRAPEZOIC INTEGRA COC See 
Se =. 1 

SUMA = 0.49 

DO 80C 1 = 23s K 

SUMA = SUMA + AREA(1) 

VCALTR = (DELTAL/2-eU) * (AREA 1) + 2.0 *SUMA + AREA(M)) 
mie Os) * {VOCAL MRO yy heey VACT 

Were ( 63900) 

PORMAI(@5xA* SHALPHAs 15XA*s 4HVACTm 15X+ GHVCALTR: 15X>+ 6GheekRie, 
WRITE (629400) ALPHA: VACTs® VCALTR: FERRTR 


*2ee4at ote 
; i 
i-~ te 
' 
- = ° 
of 
ss $2) 


-— 
. > @we af : = . * i 


= eo 7 eS = 


Tn gee 


- 4-6 Sage + 
2S ee, eee as Oe « hp “~~, te eyeetee 


=-> ——_ -— Goegteet ie ©281° & eee Sa 
°. 
LG 2:20 LLL} 
oo SEF ’ = - see 
<« “ea 
af. 
atl === aie «ciere ——- «=! 
-— — 
A Fig F = fi @ 
“is >= 
» 5 @ : 
a ES 
Sieg 88) = 
—a eee - 
= . ~ee i = 
—_— SS) eee § coe 6 





Vay 


400 


eae 


O20 


960 


970 


70 
meee (lok. oes 1SKat Flees 12k. Tf iste lt ok set oser 7) 
STOCKPILE VOLUME GCOMPUTEDSeITASIMPSONS INTEGRATIONS RUE 
KK = J - 1] 
KKK = KK = ] 


ODAREA = 060 

EVaREA = O80 

PO 910 1 = eK Ka 2 

EVAREA = EVAREA = AREA ae 

DO 920 | = 35 KKKe 2 

ODAREA = ODARCA +" AR eee 

VCALSN = (DELTAL/$@.)*¥{({AREA(]) + 4-*#EVAREA + 2.*ODAREA + eee oe 
PERReN = 100.0 * (VGAlIbahe Cl) 7 VACT 


Meese  ( 6s 950) 
FORMAT(15X» SHALPHA> 15X%s GHVACT> 15X>s OHVCALSN:»> 15X%s GHPERKS\ 7m 
WRITE (6240) ALPHA® VACT+ VCALSNs PERRSN 


STOCKPILE VOLUME COMPURNEReWITRREDBEES INTEGRATION RUE. 

Sy eae le 

DO 96C IT = 25 Ms 6 

SUMTWA = SUMTWA + 1-*AREA(T-1) + 5e*ARCA(I) + 1-*AREA( IT +1) + 


eee ARE A( I+2) + 1 o*AREReeee eee tReet 4) SS LA REA ee 


VCALWD = «2*DELTAL*SUMTWA 

PERRWD = 10U.0#®(VCALWD -— VACT)/VACT 

WRITE (69970) 

FORMAT(15Xs SHALPHA» 15X% 5 4HVACTs 15Xs 6HVCALWDs 15X» 6HPERRWL/. 
WRITE (629400) ALPHAs VACTs VCALWD» PERRWD 

eee O79 

Ce EAT T 

END 








PTZ 
- _ rr 
— - Fl ———— 


FLOW DIAGRAM FOR COMPUTER PROGRAM TO DETERMINE THE APPROXIMATE 
VOLUME OF A STOCKPILE ASSUMING THE CROSS SECTIONS TO BE ISOS- 
Geeeeom in TANGLES. 


Stake 


V/ 
READ 
XRGiee es 1M 






Mi Gioe 










READ 
L, W, ALPHA, N, VACT 






VW? 
REAL L, N 
| 


e 
RN) De ER Go . OE eo 


ALPHAR = 3.14159* ALPHA/180.0 












> . LIB Fi PO, OFFA A wl GLEN eh © RL COP, eal Baler ABP AGY FRAP DE 2 eel lable bee ental acetate, 


Vv 





C = (SIN(ALPHAR)/COS(ALPHAR))/4.0 





ETP OOS OG Pe Ry een — z 


DELTAL ="L/ (ei) 


@) 











AREA(I) = C * BASE * BASE 


SE ety ADS LO IO | 


het hee AP Chel Vente Birla SOT le Be na 6 Cerca td ede ate in Ow bee ee + Ae 





t2 











‘ 
WRITE 
AREACT) 







ee UP TT Et em a y [OT ET ea glean eg Pag TR (TEMS OA SAORI. SA, NL AES Ae Arent Vall Ts I eB Te IN tg 


VOLUME COMPUTED USING TRAPEZOIDAL RULE 


HS 
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©) 


= | 


| V7 

| 

| SUMA = SUMA + AREA(CT) 
’ . 


OR OP Pees re wr ne matt 






PERRTR = 100.0*(VCALTR - VACT)/VACT 


Wf 
WRITE 
ALPHA, VACT, VCALTR, PERRTR 


a 


PT Tag One ee 


WY 
2 








VOLUME COMPUTED USING SIMPSON'S RULE 








(es. 











EVAREA = EVAREA + AREA(T) 
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eS Om eee 8 ee He 





| = ODAREA + AREACT) 


an 


ae 


VCALSN=(DELTAL/3. )* (AREA(1)+4. *EVAREA+2. *ODAREA+AREA(M)) 


; 


N 
ne en ee eee eee f= = eco Va a 


a -— + 














ge ee i OR 





SLE EEO eG, AD AE, 


PERRSN = 100.0*(VCALSN - VACT)/VACT 


| a 


WRITE 
ALPHA, VACT, VCALSN, PERRSN 


ee ea SE SRN Rn SE ‘> _- 


¥ a ENE EN EE AI a TD A UNE 
iin a, 















VOLUME COMPUTED USING WEDDLE'S RULE 











SUMTWA = SUMTWA + 1. *AREACI-1) + 5. *AREACI) + 
1. *AREACI+1) + 6. *AREACI+2) + 1. *AREACI+3) + 
5. *AREACI+4) + 1. *AREACI+5) 










WRITE 
ALPHA, VACT, VCALWD, PERRWD 





Ww 
ay Me 


| CALL LxXas | 


a 





ES 


COMPUTER PROGRAM TO DETERMINE THE APPROXIMATE VOLUME OF A 
STOCKPILE BY DOUBLE NUMERICAL INTEGRATION UTILIZING THE 
TRAPEZOID, SIMPSON'S AND WEDDLE'S QUADRATURE RULES. 


Purpose 


The purpose of this program is to determine the approxi- 
mate volume of a stockpile by numerical integration of the 
height measurements for each cross section using the Trape- 
zoidal, Simpson's and Weddle's Rule. Then, using these re- 
sults, determine the volume of the stockpile through the use 
of the three quadrature formulas. Throughout the calculations, 
the number of height measurements and the number of stockpile 


cross sections are designated as parameters. 


LANGUAGE 


Fortran IV CIBM 7040 Computer ) 


SYMBOLIC DICTIONARY 


VARIABLE Sa 1/O*xx DESCRIPTION 
CRSECT A i Height measurements in feet 


of the stockpile. 


AREATR A O Computed area of cross sections 
using the Trapezoid Rule. 


AREASN A O Computed area of a cross section 
using Simpson's Rule. 


AREAWD A O Computed area of a cross section 
using Weddle’s Rule. 


*S - Single Variable; A = Array of Variables 
*xT = Input; O — Output 





ae 1 


DIFF 2 


Bier 3 


N PW 


NPL 


N PWX 


NPLX 


W 


VACT 
DELTAW 
BELTAL 


SUMPTS 


ODD PTS 


EVNPTS 


Yn Hn HA WN 


1/O 


1/o 
1/0 


I/O 


19 


Difference in area computed by 
the Trapezoidal Rule and Simpson's 
Raa le 


Difference in area computed by the 
Trapezoidal Rule and Weddle's 
Rute: 


Difference in area computed by 
Simpson's Rule and Weddle's 
Rule. 


Maximum number of cross section 
height measurements. 


Maximum number of cross sectional 
areas for the stockpile. 


Actual number of cross section 
height measurements for a 
specific calculation. 


Actual number of cross sectional 
areas for a specific calculation. 


Distance between the end 
boundary lines. 


Distance between the side 
boundary lines. 


Actual volume of the stockpile. 
Cross section interval size. 
Length interval size. 


A holding variable used to 
accumulate the height measure- 
ments for calculating the areas 
of cross sections with both the 
Trapezoid Rule and Weddle's Rule. 


A holding variable used to accu- 
mulate the odd numbered height 
measurements for calculating the 
areas of cross sections with 
Simpson's Rule. 


A holding variable used to accu- 
mulate the even numbered height 
measurements for calculating the 
areas of cross sections with 
Simpson's Rule. 





SUMA 


ODAREA 


EVAREA 


VCTRTR 


VCTRSN 


VCTRWD 


VCSNTR 


VCSNSN 


VCSNWD 


VCWDTR 


VCWDSN 


VCWDWD 
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A holding variable used to accu- 
mulate the areas of cross sections 
for calculating the volume of the 
stockpile utilizing. the teaee. 
zoidal Rule and Weddle'’s Rule. 


Apnholding vardable used to -aceus 
mulate the odd numbered areas of 
cross sections for calculating 
the volume of the stockpile 
utilizing Simpson@s Rule. 


A holding variable used to deers 
mulate the even numbered areas 

of cross sections for calculating 
the volume of the stockpile 

uti ligane Sampson’ s Riuiiece 


The volume of the stockpile com- 
puted by numerical integration 
of the cross section height 
measurements by the Trapezoid 
Rule and then numerically “inre— 
pratang thesresulbtine  croscecees 
tional areas Dygrne lrapezerd 
Rule to obtain the volume of the 
stockpile. 


Volume computed using Trapezoid 
Rule and Simpson's Rule. 


Volume computed using Trapezoid 
Rude anid Weddle sit le. 


Volume computed uSing Simpson's 
Rule and Trapezoid Rule. 


Volume computed using Simpson's 
Rule twice. 


Volume computed using Simpson's 
Rule andwWeddie's Rute 


Volume computed using Weddle's 
Rule and Trapezoid Rule. 


Volume computed using Weddle's 
Rule and Simpson's Rule. 


Volume computed using Weddle's 
RG sey eer 
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PLInTR S O mhiemper Cent €rrorn in Thieme. 
culated volume when the calculated 
volume is obtained by numerically 
inteprauing the cross sect nom 
height measurements with the 
Trapezoid Rule and subsequently 
numerically integraring | pieeues 
sulting cross sectional areas 
USING et he: dcapezoiden ue. 


PETRSN S O The per cent error in the wear 
culated volume utalizing Trape= 
zoid and Weddle's Rules. 


PESNTR S O The per cent error in the cal- 
culated volume utilizing Simp- 
son's and Trapezoidal Rules. 


PESNSN S O The per cent error in calculated 
volume utilizing Simpson's Rule 
EWaGe- | 

PEWDTR > O The per-cent error in calculated 


volume utilizing Weddle's and 
Trapezoidal Rules. 


PEWDSN S O The per cent error in calcwlaresa 
volume utilizing Weddle's and 
Simpson's Rules. 


PEWDWD S O The per cent error in calculated 
volume utilizing Weddle's Rule 
twice. 


PROGRAM ROUTINE 

Thas program utilizes as data the height measurements, 
at one inch intervals, of the stockpile cross sections. The 
area of each cross section is then computed utilizing all 
three quadrature formulas while varying the number of height 
measurements from the maximum to the minimum. 

The resultant areas are then numerically integrated by 
three quadrature formulas to compute the volume of the stock- 
pile. The number of cross sectional areas is varied from the 


maximum to the minimum during this phase of the computations. 
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The resulting nine combinations for calculating the 
stockpile volume are then compared to the known volume to 


obtain the per cent error in the calculated volume. 
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Does PROCRA CALCULATES! Tet VOCUMPesk AN [KRESULAR SRUCKE Ec ~ 
VOUBLE INTCGRATION USING THE TRAPEZOIDALs SIMPSONS» AND 
MEODLES RULES. 


PeeINNSITON CRSECT (305/o)s AXEATRI Vays AREASN(MO)s AREAWOaaEee 
peeeNSION DIFFI070) s’ DIF FZ rete Fr 5170) 

READ (545) News NeL 

FORMAT (213) 

WRITE (634) 

FORMAT (11) 

feeD (5326) NPWXs NPLX 

FORMAT (213) 

WRITE (636) 

meemaAl(1SXs22HN0. OF CRSECT DATA RIT Ss5Xs22HNO~e CF LENGIM DA ee 
WRITE (6097) NPWA 9 NPLX 

MeRMAT (25X%s 15s 22Xs I13/7/) 

fet 5510) (I CRSECT(1+J)s Tate Wew)>s sels Nee) 

momma T «(1 OF SaeB/7TF5D.3) 


CONVERTING DATA VALUES INTO INCHES 


Bee J = ts NPL 
DO 11 I = ls NPw 
Meoeel fist) = CRSECT(1sJ)*1l2eG 


fie) (520) Ls Ws VACT 

memrAl (2135 Fl0e2) 

WRITE (6221) 

mere! (15Xs ZOHINPLT DATA CONSTANT S77) 
WRITE (6322) 

BmOrMAr (15X5 lHL» 15X>s 1HWs 15Xs 4HVACT /7) 
Viobbeeioe23) Ls» Ws VAGH 

meer (i 5X>+ 13+ 12X>s 13s BX: FlOe2Z/7/) 
INTEGFR WwW 

DELTAW = W/(NPWX —- 1) 

DELTAL = L/(NPLX —- 1) 

WRITE (6324) DELTAWs DELTAL 

mee) (lox >s GHDELTAW = » F5e22 5X8 SADELTAL = » F5.277) 


AREA OF CROSS SECTIONS COMPUTED BY TRAPEZOIDAL» SIMPSONS 
AND WEDDLES RULES 


Mmels TRAPEZOID RULE 


—— NPW = 1] 

Bao 30e) = 1» NPL 

SeMPTS = 3-0 

BS 31 | = 25 K 

ries = SUMPTS + CRSECT(IsJ) 

Beer dR (J )={DELTAW/2-¢)*(CRSECT(1sJ)4+2e*SUMPTS+CRSECT (NPWsJ)} 


PDs SIMPSONS RULE 


Kke = NPWe- 4 
KKK = KK - 1 
Pee = 1s NPL 
ODDPTS = 0.0 





CY CxO?) 


COV CY Ar) 


errr et 
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Were t = 2+ KKse 2 

Seer = EVNP TS + CRSECT (lad) 
DO 344 M = 3 KKK 6 2 

Some oe ta. = ODDPTS + CRSEC Mee) ) 

Scene) = (DELTAWZS. PX CREECT Oise) + 4e*EVNPTS 4°25 -ODDP IS 0+ 
Pepe Ti News J) } 


PRUs WEDDLES RULE 


Pemec J =.15 NPL 
SUMPTS = Ved 
meee | = 25 NPWs 6 
eee 1S = SUMPTS + le*¥CRSECTUI—-1sJ) + SeeCRSEGIR (lea) 
fae CRSECT( 14leJ) “Mose chSeCi(iltesJ) + le*CKSECT (13s 
fee ROLCCT( 144355) + Let CRSECTUI1453 J) 
38 AREAWD(J) = « 32*SUMPTS*DELTAW 
Boe 7 | = te NPL 


DErPPY{1) = AREASN(1) - AREATRYT) 
PAPR2(i1) ="AREAWD( 13 — AREA) 
25 DIFF3(1; = AREAwWD(I) - AREASN(T) 


Pelle (6935) 

Moe FORMAT(15X» GHAREATR:e 10X2 GHAREASNs 10X22 GHAREAWDs 10X%s SHDIFFI> 
Meee SDHDIFFAm OX>s SHDIFF3/7} 
meet iE (6536) (le AREATRI1)= AwpeaSN( ls ARCAWD(1) »DIFEP tl. 
meme r ec (li}s DIFF3Z(1)s 1 = eae) 

Somer MATIC LOXs I135 Fl0e4s 6GXs FPlOse4s OXs FlOCe4s 5Xs Fleas) Gye? os 
eee Fe 4} 
WRITE (6:39) 

39 FORMAT (3X) 


COMPUTATION OF VOLUME USING TRAPEZOIDAL»s SIMPSONS AND WEDDLES 
RULES. 


Mote AREATR AND LENGTH USING TRAPEZOIDAL RULE 


N = NPL - 1 
SUMA = Gel 
DO 40 JT = ?4 N 
40 SUMA = SUMA + AREATR(1I) 
wCTRTR (DELTAL/2.0) * (AREAL RCTs 2eeSUMA + AREATR (NPE) 
EETRTR LOO eO*(VCTRER — VAG Tee, C7 


mimes AREATR AND LENGTH USING SIMPSONS RULE. 


ieee NPL - 1 


NN = N - 1 
POAREA.=.0.0 


DO 50 I = 25 Ne 2 
50 EVAREA = EVAREA + AREATR(1I) 
DO 51 1 = 33 NNe 2 
PemeOoRREA = OCAREA + AREATR(1) 
ee oN = tDELTALY Se )*(AREATR(1) + Ge*#EVARCA + 2e*ODAREA + 
LAREATRINFL) ) 
Petron = POU. I*(VCTRSN = VACT)/VACT 
Bees waAkenoN AND @MENGTH USING TRAPEZOJDAL RULE 





Las —_ -_ = - 
oe * 
8 
- 
EE — —- —_ -_—=— 
® 
~_—_—-—<- =-— =| — rrr ee -~-_-— =- 
e * 


6 e = a ¢ 
ee ee — a 
- - ll ge 
° 
— —_ 








MOA 


rayne 


yo} 


Nm’ 
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N = NPL - 1 


SUMA = 0-0 
Pons! = 2s N 
60 Geer = SUMA + AREASNIT) 
Weave = (DELTAL/26)*(AREASN( 1). + Gee SUM + AREAS ee 
PESNTR = 1L1600-60*(VCSNTR - VACT)/SVACT 


4TH» AREASN AND LENGTH USING SIMPSONS RULE 


Dime NEL eel 


NN = N - 1 
Saxe A = 0.0 
EVAREA = 0.20 
eo 70 | = 25 Ne 2 
Meme eREA = EVAREA + AREASN(1) 


DO 71 I = 39 NNo 2 
11 ODAREA = ODAREA + AREASN(T) 
VCSNSN = (DELTAL/3«¢)*(AREASN(1) + 4e¢#EVAREA + 2e*ODAREA + 
LAREASN(NPL) ) 
PESNSN = 100¢0*(VCSNSN = VACT)/VACT 


Pig. AREATR AND LENGTH USING WEODEES RULE 


SUMA = Qed 
Poe?7> | = 2s NPLs 6 

foe SUMA = SUMA +l etARGATR ( lemde) cee Se*ARGATR(1) +15e*AREATR(I+1) + 
166¢*AREATR( 142) + Le*ARCATR( 143) + 5e¢*AREATR(IT+4) + 1¢*AREATR(I+5) 
VETRWD = «3*SUMA*DELTAL 
PETRWD = 100-0*(VCTRWD - VACT)/VACT 


Silj® AREASN ANDO LENCTH UteING WEDDLES RULE 


Sov A = Quid 
MOv7S | = 23 NPL; 6 

MommouMA = SUMA + 1e*ARCASN(I~-1) + 5e*#AREASN( 1) + 1.*ARFASN(I4F1) + 
Mew ARE ASN( 142) + 1-e*#AREASN( 143) 4+ SeeAREASN(I1+4) + 1.2.*AREASNGI Sy 
VCSNWD = «3*SUMA*DELTAL 
BESNWDe = LOUeC*#(VCSNWD — VACT) /VAGT 


7THs AREAWD AND LENGTH USING TRAPEZOIDAL RULE 


Meee NPL =- i 
SUMA = 020 
DO 77 I = 25 N 
77 SUMA = SUMA + AREAWD(1) 
VCWDTR (DELTAL/2¢)*(AREAWD(1) + 2e¢*SUMA + AREAWD(NPL) ) 
PEWDTR 1O00-C*(VCWDTR -VACT )/VACT 


Batis AREAWD AND LENGTH USING SIMPSONS RULE 


N = NPL = ] 


NN = N- 1 

ODAREA = 0-0 
EVAREA = 0.0 

Bo 78 [ = 25 No 2 


Bee EVAREA = EVAREA + AREAWD(1) 
woo 1 = 345 NNs 2 











avery 
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79 ODAREA = ODAREA + ARFAWD(T) 
VCWDSN = (DELTAL/3-.)#(AREAWD(1) + 4¢*EVAREA + 26*ODAREA 
1+ AREAWD(NPL) ) 
PEWDSN = 100-O0*%(VCWDSN - VACT)/VACT 


9THs AREAWD AND LENGTH USING WEODDLES RULE 


SUMA = 0e0 
DO 80 I = 23 NPL» 6 

80 SUMA = SUMA + 1-*AREAWD( 1-1) + 5¢*AREAWD(1) + Le*AREAWD(I+1) + 
166¢*AREAWD(1+2) +1.%AREAWD(1+3) + 5e*AREAWD(1+4) + 16*AREAWD(1+5) 
VCWDWD = « 3*SUMA*DELTAL 
PEwOWD = 1002¢0*(VCWDwD — V4ACT)/VACT 


WRITE(6941) Ne oT 
See ORMAT (53)Xs LOHVACTUAL = we Faiew277 } 
WRITE (6242) 

42 FORMAT(13X=s 6GHVCTRTR»s 6X2 GHVCTRSNs 6X9 OGHVCSNTRK:s 6X2 GHVCSNSN» 6X 
le GHVCTRWDs 6X5 6HVCSNWDs 6X* GHVCWOTR: 6X2 6HVCWDS\s 6X» 6BHVCWDWL 
pee OxX// ) 

WRITE(63952) VCTRTRs VCTRSN»s VCSNTR:es VCSNSN»s VCTRWDs VCSNWD » 
IVCWDTR»e VCWDSNs VCWDWD 

Meer ORMAT(LOX>s FideZ2s 2Xs FiGePemmers PlOs2s 2Xe FlOse2s 2XS er 1 One eee 

Mel Ue2>s 2Xs Floeas ZAs FYUseZs mere 10a2s 2X77} 
Meat (6562) 

62 FORMAT (13X9 6HPETRTR»e 6X9 6HRETRSNs 6X2 6HRESNTRse 6X9 
LOGHPESNSNs 6X9 6HPETRWDs 6X9 6HRPESNWDs 6X2 GHREWDTRs 6X2 GHRPEWDSN:> 
26Xs BHPEWDWDs> 6X//) 

WRITE(6972) PETRTR»s PETRSNs PESNTR:e PESNSN» PETRWDs PESNWD> 
IPEWDTR:s PEWDSNs PEWDWD 

Mee ORMAT(1L3X>s F6e2s 6X9 Fels OXeuF6e25 BXs Foe2s Gis F6eZ2s Gx 

PeO*e2 > 6X F6e2s 6X Foes: OX Peez77 } 


CALL EXIT 
END 











len 
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FLOW DIAGRAM POR COMPUTER PROGRAM TO DETERMINE THE APPROXIMATE 
VOLUME OF A STOCKPILE BY DOUBLE NUMERICAL INTEGRATION UTILIZING 
THE TRAPEZOIDAL, SIMPSON'S, AND WEDDLE'S QUADRATURE RULES. 





CRSEGT (ieee. 0 


CRSECT (€25-.)9 = on 


Se 


7 


READ 
NPW, NPL 


M 
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NP2 = NPW - 1 





READ 
CCCRSECT (C1,J), 1 =e — leer) 


t 






CONVERTING DATA VALUE INTO INCHES 





CRSECT (1,J) = CRSECmE@=n) = 12.0 





eo 












READ 
L,W,VACT 


INTEGER W 


DELTAW = W/(NPW - 1) 
io 


DELTAL = L/(NPL - 1) 


WRITE 
DELTAW, DELTAL 


AREA OF CROSS SECTIONS COMPUTED BY 


TRAPEZOIDAL, SIMPSON'S, AND WEDDLE’S RULES 
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FIRST: ‘TRAPEZOIDAL RULE 





| 
SUMPTS = SUMPTS + CRSECT (I, J) 
| 
: 


. 
CD pb I = I +1 
> 


ee - 





le 


AREATR( J) = CDELTAW/2-) eae CRS ECL as 
2. * SUMPTS + CRSECT (NPW, J)) 





vi 








SECOND: SIMPSON'S RULE 


EVNPTIS = EVNPTS@eecRSect (kK. J) 


V/ 
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AREASN(J) = (DELTAW/3.) * CCRSECT(1,J) + 
4. * EVNPTS + 2. * ODDPTS + CRSECT(NPW, J) 


ior 
re = tle 
“™ 
7 


THIRD 
WEDDLE SS] RUBE 
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SUMPTS = SUMPTS + 1.*GRSESEGI=9 J 5] Chol Glh@l nome 
1. *CRSECT(I+1,J) + 6.*CRSECT(I+2,J) + 1.*CRSECT(I+3,J 


+ 5.*CRSECT(I+4,J) + 1.*CRSECT(I+5,J) 


i Ne ee ee 






madi 





AREAWD(J) = .3 * SUMPTS * DELTAW 





94 


vine ~ AREATR(I) 


DIFF 1 (I) 





OS 
» 


AREAWD(I) ~ AREATR(I) | 
| DIFF. DIFF 3 (I) = AREAWD (I) ~ AREASN(I) | 


| DIFF 2 (I) 








nt aR TE 5 EE SA PRR 





| 
: ? Leelee i} 


™~ 
? 
WRITE 


(I, AREATR(I), AREASN(I), AREAWD(1I), DIFF 1 (I), 


DIFF 2 (1), DIFP 3 (1), I = 1,NPL) 









COMPUTATION OF VOLUME USING 
TRAPEZOIDAL, SIMPSONS Nee DEE «s 
RULES 






FIRST: AREATR AND LENGTH TRAD 











oS 







SUMA = SUMA + AREATR (TI) 





VCTRTR = (DELTAL/2.0) * (AREATR(1) + 2.*SUMA + 
AREATR(NPL) ) 


' 
e 
| mg ee RST Aap a Nr A, SRA eT A - 


ed 


100.0 * (VCTRTR-VACT) /VACT 


SECOND: AREATR AND LENGTH USING 
SIMPSON'S RULE 
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EVAREA = EVAREA + AREATR(I) 





97 








le 


\y 
| ODAREA = ODAREA + AREATRCI) | 








VCIRSN = (DELTAL/3.)*(AREATR(1) + 4.*EVAREA + 


2.*ODAREA + AREATRC(NPL) ) 


re 


PETRSN = 100.0((VCTRSN-VACT)/VACT 





THIRD: AREASN AND LENGTH TRAP 








ve 





SUMA = SUMA + AREASN(CTI) 


PESNTR = 100.0*(VCSNTR-VACT)/VACT 
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ODAREA = 0.0 





EVAREA = EVAREA AREASN (TIT) 


- 
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i 
VCSNSN = (DELTAL/3.)*CAREASN(1) + 4.*EVAREA + 2.*ODAREA + 





AREASNC(NPL) 


PESNSN = 100.0*(VCSNSN = VACT)/VACT 







FIFTH: AREATR AND LENGTH USING WEDDLE'S RULE 


New 
ee ae ere ee eee 


SUMA = SUMA + 1.*AREATR(I-1) + 5.*AREATR(I) + 
1.*AREATR(I+1) + 6.*AREATR(CI+2) + 1.*AREATR(CI+3) + 
5.*AREATR(I+4) + 1.*AREATRCI+5) 





OE a ee eee (Oe a a mn ER ee ee 


VCIRWD "Seo CoUMAS ae ELIAL 





— 


PETRWD = 100.00 * (VCTRWD-VACT )/VACT 


2 


eg ae 
- oo ss 


AREASN AND LENGTH USING WEDDLE'S RULE 






SPEXTH: 





a A ee Emr eS SAN A me 





Or 


ae. pe a aes ee Ce! OE ee Re 


SUMA = SUMA + 1.*AREASN(I=-1) + 5.*AREASNCI) + 
L.*AREASNCI+1) + 6.*AREASN(CI+2) + 1.*AREASNCI+3) + 
5.*AREASN(I+4) + 1.*AREASN(I+5) 


V/ 








7 





Oe: 


i 


VCSWD = .3*SUMA * DELTAL 


8 tee ee 





PESNWD = 100.00 * (VCSNWD - VACT)/VACT 






SEVENTH: AREAWD AND USING TRAP RULE 


TS Lem ha A ET 
spell aa 











POs 











VCWDTR = (DELTAL/2.) * CAREAWD(1) + 2.* SUMA + AREAWD(NPL 


EIGHTH: AREAWD AND LENGTH USING SIMPSON’'S RULE 


EVAR Ee 


I 
SS 
© 
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VY 


EVAREA = EVAREA + AREAWD (1) 





ODAREA = ODAREA + AREAWD (1) 


VCWDSN = (DELTAL/3.) * C(AREAWD(1) + 4. * EVAREA + 


2. * ODAREA + AREAWD(CNPL) ) 


PEWDSN = 100.00 * (VCWDSN-VACT )/VACT 
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NINTH: AREAWD AND LENGTH USING WEDDLE'S RULE 







ER a ere ah IRI et ee seo ——= 


SUMA = SUMA + 1, * AREAWD(I-1) + 5. * AREAWD(I) + 
1. * AREAWD(I+1) + 6. * AREAWD(I+2) + 1. * AREAWD(I+3) + 
5. * AREAWD(I+4) + 1. * AREAWD(I+5) 











VCWDWD = 23) * SUM ape Dee 


PEWDWD = 100.0 * (VCWDWD - VACT)/VACT 
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WRITE 
VCTRTR, VCTRSN, VCSNTR, VCSNSN, VCTRWD, VCSNWD, VCWDTR, 
VCWDSN, VCWDWD 








WRITE 
PETRTR, PETRSN, PESNTR, PESNSN, PETRWD, PESNWD, PEWDTR, 
PEWDSN, PEWDWD 


CALE Eid 





LOF 


Undisturbed Stockpile Depth Measurements 
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COMPUTER PROGRAM TO DETERMINE THE APPROXIMATE VOLUME OF A 
BANKED STOCKPILE ASSUMENG THE STOCKPILE CROSS SECTIONS ARE 
IbCserhLES TRIANGLES. 
lEerpose 

The purpose of this program #8 to @ecermine the appre. =. 
maut® VOolune of &@& #@TOcCKHL1G whobe BAGdES are in EattAct wath Tex 
Baining walis. A secondary purpose 1S to Cetermane tae anete 


of repose of the material in the stockpile. 


Language 
FORTRAN IV CIBM 7040 COMPUTER) 


Be mpolic Dictionary 








VARIABLE S/A* I /O** DESCRIPTION 
PW S I . Distance between apex of the 
Stockpile and the réevatnine ae 
WE A -= Width of the small triangular 
section Cut off by the retaistwne 
wall. 
XR A I Distance the base of the pile is 


from right boundary line. 


XL A I Distance the base of the pile is 
from the iert boundary line. 


M S I/O Number of measurement points. 
ie S I/O Distance between the end boundary 
JI neSe 
W S 1/O Distance between the side bound- 
ary” Paes 
ALPHA S 1/O Inclination angle of the stock- 
pile - same as apparent angle 


of repose. 


*S - Single variabie; A - Array of variables 


**T - Input; O - Output 





Lie 


AL PHAR -- Alpha expressed in radians. 


N I Number of measurement stations. 
_8 

VACT I/O Actual’ volume of the stockpile. 

PELIAL Leneth ism8erval sige. 


BASE = Width of stockpiie at the base. 


- WM MM HM BH 
S 


AREA O Area of a cross section. 

SUMA S -- A hoiding variabie used to 
accumulate the areas of cross 
sections for calculating the 
yotume Dyeitee Trapezoidal irer 


VCALTR = O The caiculated volume of the 
STOCGKDLLe UitiliZzine the Araces 
ZO Cae e. 


PERRTR S O Per cent error Of the caleul ce 
volume Ug@ne the Trapezoidai Ries 


table‘used to accu- 
eas or tne odd num 
le cross sections 
for calculatine oe stecaiee 
volume using Simpson's Rule: 


ODAREA S -- A holdisg va 
mulate the a4 
pered stocks 


BVAREA S -- A holding variable used Toracen= 
mulate the areas of the even 
numbered cross sections for cal- 
CUlaAEINee Ae STOCK 21Le volume 
USINegeiApson s Rule. 


VCALSN > O Lhe calculféed stockpile volume 
using Simpson's Rule. 


PERRSN S O Per cent error of the calculated 
volume using Simpson's Rule. 


SUMTWA S -— A hoiding variable used to accu- 
mulate the weighted areas of 
cross sections for calculating 
the steckpiie volume using Weddle's 
Ruves 


VCALWD 5 O The caLith seared stockpile volume 
using Weddiets Rule. 


PERRWD S O Percent error of the calculated 
Volume WwSans "“Wsddie*sS Rule: 





Pe 


Program Routine 


ithe program utilizes as datapoints the eGictanees 
which the stockpile cross sections are from the side boundary 
lines. The area of each cross section is then computed assum- 
ang that at 1s an isosceles triangle. The program checks to 
see if the stockpile is in contact with the retaining wall. 

If it is, then the program calculates the net area of the 
cross section by suber cede ome area of the cut-off tri- 
angular section from the gross area of the standard isosceles 
eraangiles. 

These areas are then antegrarea msine the trapezoracaus 
Simpson's and Weddle'*s quadrature formulas for a set inter- 
val size. The error in the calculated volume is then conm- 
muLec . 

The angle of repose is determined when the calculated 
volume error is a zero value based upon varying the parameter, 


BoC PHA . 
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THIS PROGRAM DETERMINES THE APPROXIMATE VOLUME OF A STOCKPILE 
WHICH ITS BANKED AGAINST A RETAINING WALL AND FURTHER 
me OUMING THAT THE CROSS SECTIONS ARE TSOSCELES TRIANGLES. 


DIMENSION ~XL(100)s5 XR(100)+- AREA(100)s BASEO2(109)5 WF(100) 
READ(5326) PW 

mOoRMAT TFS. 1) 

RED (5.5) M 

meRimraAT (12) 

WRITE ( 6396) 

FORMAT (15X»s 20HM = NOe OF DATA PTSe//) 

WRITE (697) M 

PORMAT (19X%s I2/7/) 

Rene 20) (AL( I) » XR( I) »s 1] =e my 

mOmRMAT (2F100e2) 

WRITE (6611) 

FORMAT (15Xs I7HINPUT DATA POINTS//) 

WRIWE (6513) 

FORMAT (15X» 3HNOes 15Xs 2HXLs 15Xs 2HXR//) 

TER (6512) ( ITs XL(I)s XR(I) » J] = le M) 

poerrm ( 15X» 125 10Xs Fl0e2; Bree F1O062//) 

Pome (5510) Ls We ALPHAs Ns VACT 

POMMAT (2F10«2sF66.42 F602 FI0.7) 

WRITE (6914) 

FORMAT (15Xs 2OHINPUT DATA CONSTANTS//) 

WRITE (6916) 

FORMAT(15X»s 1HL»>» 10X»s IHWs 10X+ SHALPHA>s 10X%s IHN» 10X%+ 4HVACT//} 
wer TE Noel?) Ls Ws ALPHA. Nee VaAcT 

SGMAT( 1 3X>. Pe2s6Xs FSe2> SMeete See2ne8 X eee 569 Xs ee 2 7 7) 


pooem., CROSS SECTIONS ARE ISOSCEUES TRIANGLES FOR COMPUTING ARreeee 


REAL LoN 

ALPHAR = 3214159 * ALPHA / 180.0 

A = SIN(CALPHAR)/COS(ALPHAR) 

mer AlL ="L/UN-1.-0) 

min] (Be 6 s 198) 

mere (15%) ITHINTERVAL = DELTALY7) 
WRITE (6519) DEAL 


MORMAT (22X59 F5e2//) 

Jo= N 

DO 44 IT = 1lsJ 

Mrmr ( 1)e«LEeO0e-0) GO TO 33 

PpaoeO7(l) = (W—-XR(I)-XL(1))72.-.0 

pmmetl) = (BASEO2 (1 )#*#2 742A 

sO, TO 44 

BASEO2(1) = (W-2.0#XL(1I)) /2.0 

WF(I) = BASEO2(1) -—Pw 

AREA( IT) = (BASEO2(1)**2)%A -— (WF(I)**2)*A*O.5 
J = | 

CONTINUE 

Peat!) = 0.0 

AREA(J) = 0.20 

Weehe (63500) 

merrmeat (15X>+ 22HAREA OF CROSS SECTIONS//) 
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Meee (63600) (is AREAVIT)s 1 = Le 'M) 
BeRMATL(LISXs Il2s FlOea) 
STOGKPILE VOCUME COMPUTED” WITH THE TRAPEZOID INTEGRATION RU ila. 


aoe «(i 

SUMA = 0.20 

Poeo00 1 = 25 

SUMA = SUMA + AREA(I) 

VCALTR = (DELTAL/2-0) * (AREA(1) + 2¢€0 *SUMA + AREA(M)) 
PERRTR = 100«e0 * (VCALTR = VACT) / VACT 


WRITE ( 69900) 

FORMAT (15Xs SHALPHAs 15Xs 4HVACTs 15Xs 6HVCALTRs 15Xs 6HPERRTR//) 
HeItie (6s4C0C) ALPHAs VACTs VCALTRs PERRTR 

MORMAT (15Xs FSe2s 15X 6 F fle 2uuio Xs Fes Coe oe 


STOCKPILE VOLUME COMPUTED WITHSIMPSONS INTEGRATION RULEe 


r= J = /1 

KKK = KK = j 

ODAREA = 0.0 

EVAREA = 0.0 

DO 910 [| = 25s KK» 2 

PyVAREA = EVAREA + AREA {17} 

DO 920 I = 35 KKKe 2 

ODAREA = ODAREA + AREAs(1) 

VCALSN = (DELTAL/34)*{( AREA] ) +9" EVAREA + 2-*ODAREA + AREA(M)) 
PERRSN = 10Ue0 * (VCALSN - VACT) / VACT 


WRITE ( 69s 950) 
FORMAT(15Xs SHALPHAs 15X»s 4HVACTs 15Xs 6HVCALSNs 15X» 6GHPERRSN//) 
WRITE (6540C) ALPHA s VACT » VCALSNs PERRSN 


SITOCKPIEE VOCUME™ COMPUTED Wl lm DOLES UINTEGRA] Tey Ue 
SUMTWA = 0-0 

BO 950 [1 = 25 "Ms 6 

SUMTWA = SUMTWA + 1e*ARCA(T-1) + 5e¢*AREA(I) + 1Le*AREA(I+1) + 


PSe*AREA(I+2) 4 1e*AREA (1+ 38 i ORR ee +) Pee eREA C1 +37) 


VCALWD = .«3*DELTAL*SUMTWA 

PERRWD = 100¢O*(VCALWD - VACT)/VACT 

Meite (69970) 

FORMAT(15Xs BHALPHAs 15Xs SHVACTs 15Xs GHVCALWDs 15Xs 6HPERRWD//) 
WRITE (63400) ALPHA» VACTs VCALWDs PERRWD 

Bo 70 9 

SAblL EXIT 

END 
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PLOW DIAGRAM FOR COMPUTER PROGRAM TO DETERMINE THE APPROXIMATE 
VOLUME OF A BANKED STOCKPILE ASSUMING THE CROSS SECTIONS TO BE 
MSOSCELES TRIANGLES. 







XLOTSaexRGio, 
I=1,M 


READ 
L, W, ALPHA, N, VACT 














ets 


\ 


DELTAL = L/(N = 1.0) 


+ 
y 





WRITE 
— 


eo 










| 


| BASE 02(I) = (W-XR(I)-XL(I))/2.0 | 





7 
fA 


Li6 





WRITE 
AREA(I) 


VOLUME COMPUTED USING TRAPEZOIDAL RULE 








Le 






SUMA = SUMA + AREACT) 







PERRTR = 100.0*(VCALTR - VACT)/VACT 


WRITE 
ALPHA, VACT, VCALTR, PERRTR 








ies 






ODAREA 






BVARLE 





EVAREA = EVAREA + AREA(CI) 
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oo 


ODAREA = ODAREA + AREA(CT) 


Fe ee etn ne 





VCALSN=(DELTAL/3. )* (AREA(1)+4.*EVAREA+2.*ODAREA+AREA(M) ) 


RR, CL Oy TN TTL BPD FCC aC Ei lt Sta ITS LOE AAO Ltn site OO eat Pas 








WwW 





PERRSN = 100.0*(VCALSN - VACT)/VACT 






WRITE 
ALPHA, VACT, VCALSN, PERRSN 
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SUMTWA = SUMTWA + 1. *AREACI-1) + 5. *AREACI) + 
1. *AREACI+1) + 6. *AREACI+2) + 1. *AREACI+3) + 
5. *AREACI+4) + 1. *AREACI+5) 


RY 
> 











-3 * DELTAL * SUMTWA 


WRITE 
ALPHA, VACT, VCALWD, PERRWD 





APPENDIX B 
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FIGURE | 
CALCULATION OF STOCKPILE VOLUME FROM BASE MEASUREMENTS 


Boundary Width | Anrle of Renose C = Tone /4,0 | Interval] Sise,J} Boundary Leneth 


sBL/(# of Cross 

















Meco. (1) | Col, (2) | oi, (3) [ Col. (4) [cod (5) (cote) cot. (7) C0 (8) 
5 


Bre Width of Cross soct. ODD Cusrli) [CeaVaN CxtucT 

















Section Pile Base Area x 2 x & 
Gol (6) xc Col (7) ODD | Col (7) evs: 
Ft4 4 











TOTAL 


HTL 


ESTIMATE USING TRAPEZOID RULE ESTIMATE USING SIMPSON'S RULE 


Volume = Interval x Total Col (7)= Volume = Interval x Totsi Col (8)+ Total Col (9)= 


Tonnage = Volume x Bulk Density x (_1.), 
(2000 ) Tonnege = Volume x Bulk Density x (_ 1.) ¢s 
(2000 ) 
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FIGURE 2 


CALCULATION OF STOCKPILE VOLUME FROM BASH MEASUREMENTS 


Interval Stxe Is |Moundary Length 
TA A eee0 Gere BL/(# of Uros3 hI = 2 N ite 


sO, 3 Col.(7i.. (“71 eel some ool een 
Cross ; i ‘ Cross Sect. Cou CASl3T) | (eve.. onoe 
Sect ion Pile B- PLlé > mt Arca x2 yh 
Col({6) x,¢ Cal(7) coo Col(7) EVEN 
x 2 


Ft. x 





ESTIMATE USING TRAPEZOID RULE {ATE USING SIMPSON'S RU. 


Volume = Interval x Tots] Col(7) = (1.0) (854.38) = 854.38 CU.FT Volume s_Interya) x(Total Col (8)+Total Col (9)) 
3 
s (1) (852.58 + 1712.30) = 854.98 CULFT. 


te es 


Tonnege = Volume x Bulk Density z (eae) Tonnare = Volume x Bulk Density x (_ 1 _) 
(2000 ) (7000 ) 
# (854 .38)(47.19)(__1_) = 20.16 TONS © (854.98)(47.19)(_1 _) « 20.16 TONS | 
2000 




















